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Abstract 

The technique of Hamiltonian flow equations is applied to the canonical Hamil- 
tonian of quantum electrodynamics in the front form and 3+1 dimensions. 

The aim is to generate a bound state equation in a quantum field theory, partic- 
ularly to derive an effective Hamiltonian which is practically solvable in Fock-spaces 
with reduced particle number. The effective Hamiltonian, obtained as a solution of 
flow eqautions to the second order, is solved numerically for positronium spectrum. 
The impact of different similarity functions is explicitly studied. 

The approach discussed can ultimately be used to address to the same problem 
for quantum chromo dynamics. 



Preprint: MPI-H-V33-1998 
PACS-numbers are: 

ll.10.Ef Lagrangian and Hamiltonian approach 
11.15.Tk Other non-perturbative techniques 

12. 38. Aw General properties of QCD (dynamics, confinement, etc.) 
12.38.Lg Other non-perturbative calculations 

*Present address: Department of Physics, North Carolina State University, Raleigh, 27695-8202 NC 
USA 

^on leave from Inst, for Theor. Physics, Eotvos University, Pazmany P.s. 1/A, Budapest, H-1117 



1 



CONTENTS 

Contents 



I Introduction 



2 Flow equations and bound state problem 



3 Canonical QED Hamiltonian on the light-front 

4 Flow equations applied to QED 



4.1 The exchange part . 

4.2 The annihilation part 



4.3 Integrating the flow equations for the exchange 



4.4 Integrating the flow equations for the annihilation 



4~5 The effective ee-interaction 



4TB Dependence on the cut-off function 



5 Discussion and interpretation 



6 Numerical solution for positronium spectrum 



6.1 Formulation of the problem . 
672" Discussion of numerical results 



7 Summary and conclusions 



A Flow equations and Feynman amplitudes 
B Denning different cut-offs 



C The matrix elements in the exchange channel 






C.l The helicity table 




C.2 The helicity table for arbitrary J z . 





D The matrix elements in the annihilation channel 



D.l The helicity table 



D.2 The helicity table for | J z \ < 1 



1 INTRODUCTION 



3 



1 Introduction 

Over the past twenty years two fundamentally different pictures of hadrons have devel- 
oped. One, the constituent quark model is closely related to experimental observation 
and phenomenology. The other, quantum chromodynamics (QCD) is based on a covari- 
ant non-abelian quantum field theory. Disregarding lattice gauge calculations, one has 
several reasons why the front form of Hamiltonian dynamics || , as reviewed recently 
in JJ, is one of the very few candidates for reconciling the two approaches. Particularly 
the simple vacuum and the simple boost properties confront with the complicated vac- 
uum and the complicated boosts in the conventional Hamiltonian theory. Wilson and 
collaborators [f|, || have proposed a scheme in which one presumes a potential for the 
bound-states and handles the relativistic effects by structures imposed by the needs of 
renormalization. The available numerical examples |j, |7| however violate admittedly some 
symmetries of the Lagrangian and it is not clear how to restore them systematically. 

There are two major problems when one addresses to solve a Hamiltonian bound state 
equation 

H\iJ>) = E\iP) (1) 

in a covariant relativistic field theory. First, the canonical field-theoretical Hamiltonian 
H contains states (fields) with arbitrarily large energies. Second, the number of particles 
in a field theory is unlimited and H contains the impact of arbitrarily many particles. An 
eigenfunction for example a meson wave function, has contributions from arbitrarily 
many Fock-space sectors \ip) = ^ Q q\qq) + <£qqg\<i<ig) + — Therefore, in general, the 
Hamiltonian operator H can be understood as a matrix with infinite dimensions both 
with respect to 'energy' and with respect to 'particle number'. The method displayed 
below, the 'Hamiltonian flow equations' of Wegner ||, cope with either of them. 

It is a subject of its own and not completely trivial to write down a suitable Hamilto- 
nian (operator) H 0. In the sequel we shall use the canonical (light-cone) Hamiltonian of 
gauge theory in the light-cone gauge (A + = 0). In the light-front Hamiltonian approach 
one faces then two classes of problems that are related with each other: the problems 
associated with the light-front formulation and the problems in formulating an effective 
Hamiltonian theory and the difficulties in the renormalization program. In order not to 
ponder all these problems with the problems of confinement and chirality, one disregards 
in this work QCD and restricts to QED as a model case. Unlike in other work [7], ||], [K| on 



QED and other models we apply in this work the method of flow equations || [12], [T3|, |H 
with the objective to derive from the Hamiltonian in the front form a well-founded effec- 
tive low-energy Hamiltonian which acts in the space of a few particles and which can be 
solved explicitly for bound states. 

Since we aim at a pedagogical presentation we sketch shortly the ingredients of front- 
form QED in Sec. |^. In Sec. § the general aspects of the Hamiltonian flow equations 
are collected in such a form that the application of the flow equations to QED in Sec. |] 
becomes more transparent. It is here where the bulk of the present work is displayed in a 
formal way. The implications are discussed in Sec. ||], and the numerical calculations are 
given in Sec. |6|. The work continues to be actively pursued, and this is only the first in a 
series of papers. 
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2 Flow equations and bound state problem 

In this work, we shall focus on the flow equations for Hamiltonians formulated first by 
Wegner || aiming at the construction of an effective bound state Hamiltonian for field 
theories. The general aspects of the method do not depend on the nature of the Hamilto- 
nian and before plunging into the paraphernalia of the field theoretical details it is useful 
to outline the general ideas in a manner slightly different from the original formulation 

!• 

It is always possible to divide the complete Fock space (with its different particle num- 
ber sectors) into two arbitrary subspaces, called the P- and the Q-space. The Hamiltonian 
matrix of Eq. (|l]) has then the form 

R= (PHP PHQ\ 

[ QHP QHQ ) ' {Z) 

with P and Q = 1 — P being projection operators. Suppose now we are unable to solve 
the eigenvalue equation for the whole matrix, because of, say, computer limitations. The 
method of flow equations allows then to unitarily transform the Hamiltonian matrix into 
a block-diagonal form 

_ / PH cS P \ 
HcS ~ [ QH cS Q ) (3) 

by analytical procedures. The two blocks are then decoupled, and the eigenvalue problem 
with an effective Hamiltonian H e g can be solved and diagonalized separately for either of 
the two spaces, which technically might (or might not) be easier than the solution of the 
full problem. 

Since one can choose the number of particles in the P-space one reduces in this way in 
general the many-body problem to a few particle bound state problem at the expense of 
finding a more complicated effective Hamiltonian operating in a limited particle number 
space. This general idea is similar to the procedure of Tamm and Dancoff [TT5], IB"] where 



an effective interaction in a few particle sector was obtained by eliminating the 'virtual 
scatterings' to the higher Fock-space sectors in the Q-space, see also f|]. 

The method of flow equations || [12], [13], [14] works with a unitary transform which is 



governed by a continuous parameter /. The unitarily transformed Hamiltonian is then a 
function of this parameter, i. e. 

^ = W),H(l))- (4) 
The generator of the transformation is subject to some choice but taken here as J8| 

V (l) = [H d (l),H(l)}. (5) 
The Hamiltonian is separated conveniently into a block- diagonal part 

HS) = ( PH j 3 l)P QH ° {l)Q J , (6) 



2 FLOW EQUATIONS AND BOUND STATE PROBLEM 



5 



and into the rest 

m-HAD = ( QH \ l)P PH « )Q ), m 

which is purely off-diagonal. It is precisely this rest which due to PHQ 'changes the 
particle number'. In many cases of practical interest one can interpret this rest as a 
'residual interaction'. If the rest vanishes, or if it is exponentially small, one has solved 
the most important part of the problem. In the sequel we convene that the flow parameter 
changes from I — to I — > oo, corresponding to a change from the initial canonical 
Hamiltonian to block-diagonal effective Hamiltonian. According to Eq.(^) the generator 
is always off-diagonal. 

~ { Qr)(l)P ) ■ (8) 

The flow equations Eq.@ for the diagonal and the rest sectors can then be disentangled 
into 

-tjPHP = Pr]QHP- PHQr]P, (9) 

^PHQ = PriQHQ - PHPrjQ , (10) 
dl 

and the trivial identity 

Pr]Q = PHPHQ - PHQ HQ . (11) 

For QHQ, QHP, and QrjP one proceeds correspondingly. Since TrPHQHP is restricted 
from above, Wegner's choice for the generator Eq.(|5|) results in a monotonously decreasing 
measure for the off-diagonal particle-number changing interaction 

^-TrPHQHP = Tr( Pr]Q(QHQHP - QHPHP)) 

+ Tr[(PHPHQ- PHQHQ)Q V P) 

= 2Tr (PrjQrjP) < . (12) 

For the generator this implies that rj(l) ^0 in the limit / — > oo and that the block-diagonal 
part of the Hamiltonian commutes with the Hamiltonian itself. Thus, the effective Hamil- 
tonian H e s becomes 'more and more block- diagonal' with increasing flow parameter. 

In general the solution of these equations will become quite involved. One reason is, 
that the equations are nonlinear. Another is, that starting with a two-particle interaction 
one generates due to the commutators three-particle, four-particle etc. interactions. It 
is however possible to solve the equations in certain limits or approximations. A limit in 
which the equations can be solved exactly to a large extend is the n — > oo limit of an 
n-orbital model 0. In this limit the equations for the two-particle interaction are closed, 
that is generated three-particle interactions do not couple back to the flow equations 
for the leading two-particle interaction. For realistic systems, which normally do not 
obey such a limit one can truncate the equations, which turned out to give very good 
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results for the Anderson impurity model [18] and the spin-Boson model |19[]. Another 



approach is to perform a perturbation expansion in some coupling. This has been applied 



to the elimination of the electron-phonon coupling [|TI|. Similarly to [13] we will use this 
approach for the positronium on the light-cone here. 

Before we enter this calculation it seems appropriate to explain, how this procedure 
works and to compare it to the similarity renormalization by Glazek and Wilson [[|. 
Suppose we would know approximately the eigenstates of the sector Hamiltonians PH(l)P 
and QH{l)Q and their eigenvalues E p {l) and E q {l). The indices p and q run over all states 
in the P- and Q-space, respectively. Suppose further, that this basis is /-independent. This 
means in other words, that we assume, the off-diagonal matrix elements h pp / and h qq > of 
PHP and QHQ are supposed to be small. Then in evaluating the commutators in Eqs. 
( [TO]) and (]TT|) we neglect the small off-diagonal matrix elements h pp i and h qq i and take 
into account only the diagonal matrix elements E p and E q . In Eq. (|9[) however we keep 
all matrix elements on the right-hand side. Then Eqs.([9|)-([TT|) yield 

dh r 



u pp , 



dl 

dh. 



dl 

Vpq 

Substituting r\ pq yields 



Yl {Vpqhqp' ~ h pqVqp>) , (13) 

q 

-(E p -E q ) Vpq , (14) 

(E p -E q )h pq . (15) 



dh^l) _ ^( dh pq {l) 1 1 dh qp ,{l) \ 

dl ~ V\ dl E p (l)-E q (l) qp ' [)+pg[) E p ,(l)-E q (l) dl )' 1 DJ 

where 

^ = -(E p (i)-E g (i))\ q (i). (i7) 

The analogous equation for h qq > is obtained by interchanging p — > q, p' — > q' . 
For the off-diagonal rest part one gets 

MO = MO) ex p( - f dl'{E p (l') - E q {l')f) . (18) 

The /-dependence in this equation can become important. If in the limit / — > oo the 
difference E p — E q vanishes, then the /-dependence can be quite crucial. If was first 
observed in the spin-Boson model fl9|] and later also for the electron-phonon coupling 
r^[, that a self-consistent solution yields a decay E p (l) — E q (l) oc 1/Vl in the case of 



asymptotic degeneracy, so that the corresponding off-diagonal matrix element h pq decays 
algebraically to zero. This procedure, however, goes beyond perturbation theory. 

One has thus reached the goal: As the flow parameter tends to infinity the rest sector 
PHQ tends to zero and is eliminated. Simultaneously, this elimination gives rise to an 
effective Hamiltonian which has the block-diagonal structure PH(oo)P and QH(oo)Q. 
The block PH(oo)P is defined by 

V(oo)-V(o)-y o dl Ep (l)-E q (l) + E p ,{l)-E q {l) dl )■ (19) 
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The first term represents the initial interaction in P-space and the second term originates 
from the elimination of the off-diagonal rest sectors. 

Here we have given a rough idea on how the flow equations work. We have to consider 
how they work in perturbation theory and we have to decide on the blocks defined by the 
projectors P and Q. Obviously the P-space should contain the states with one electron, 
one positron, and zero photons. The rest may be covered by the Q-space. Since the 
explicit calculations are done in terms of creation and annihilation operators, it is easier, 
to introduce not two but infinitely many blocks. Each block contains all states with a fixed 
number of electrons, a fixed number of positrons, and a fixed number of photons. The 
above equations can be easily generalized to this case. It is not necessary to write down 
the blocks explicitly, since the expressions in terms of creation and annihilation operators 
show explicitly, whether the number of particles is conserved or not and thus, which 
terms contribute to the diagonal and which to the off-diagonal part of the Hamiltonian. 
In quantum electrodynamics the small coupling is the charge g. The leading contribution 
is the kinetic energy of order g°. (The interactions are explicitly given in the next section). 
It yields the leading diagonal matrix elements E p , E q and no off-diagonal contributions. 
The vertex interaction, which describes emission and absorption of photons is of order 
g and purely off-diagonal, since it changes the number of photons. The instantaneous 
interaction is of order g 2 and contains both particle-number conserving and particle- 
number violating contributions. In zeroth order in g we have only the kinetic energy. 
This zeroth order contribution does not change. To first order in g we have Eqs. (|I7|JI8D , 
where the /-dependence of E has to be neglected. Then one enters Eqs. ( P^B|JH3| ), neglects 
the /-dependence of E on the right hand-side and obtains h pp > to second order in g 2 . 
The result is the effective interaction between electrons and positrons. It also includes a 
change of the one-particle energies, which become /-dependent in this order. 

Equation (|18D may be written 



We discuss now the choice of a more general 'similarity function' /. Such a more general 
function was first by Glazek and Wilson Q. 

Correspondingly, the generator of the transformation is written as 




(20) 



with 




(21) 




(22) 



An example of a different choice for / is given by 



VpqV) 



sign(P p - E q )h, 



(23) 



which yields an exponential decay of / 




(24) 
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This similarity function is good, if the sign of the difference E p — E q does not depend on 
the momenta of the interacting particles. This is the case for the absorption and emission 
of the photons in the light-cone frame and can thus be used here. Other possibilities are 
a sharp cut-off, if the energy difference is larger than a given energy of I. Glazek and 
Wilson used a continuous elimination if the energy difference lies between two energies 
which decrease with I. 

Such similarity functions can be used in two cases: 

(i) In the first case one does not introduce blocks, but aims to diagonalize single states. 
Then starting from plane waves the two-particle interaction becomes negligible, which 
may prevent the procedure from diagonalization. This happened in a first attempt in || , 
and Jones, Perry, Glazek J7| could perform the elimination of the off-diagonal interaction 
only down to energy differences of order Rydberg. This problem which shows up in the 
continuum, might be overcome, if one can introduce a discretization. 

(ii) For block-diagonalization they can be used to the order of perturbation theory as 
discussed here. If one goes beyond this order, then it is not obvious how to use a general 
similarity function. Despite the fact, that there is a lot of freedom to choose rj, one has 
to make sure, that the off-diagonal matrix elements really decay, as shown in Eq. (]12[) for 
the choice (||). 

The most important properties of the similarity function f(z) are 

/(0) = 1, 

/(2-00) = 0. (25) 

Its functional dependence on z is less important. As to be shown below in Sec. || and 
in Sec. ^] by way of example, different choices for the similarity function have almost no 
impact on physical observables like the spectrum or the wave functions of bound states. 
This leaves us finally with 

^efw = = V(°) - M°) V(°) 

* (e^wJ?^^ + .(26) 
in order g 2 . 

For p = p' this equation may contain ultra-violet divergences as I goes to 0. Remember 
that the elimination of the rest sector in Eq.fll8|) by means of flow equations is performed 



not in one step, as in the Tamm-Dancoff approach [[15] , [Tgl , but rather sequentially for 
different energy differences, i.e. 

A= i=<| £p - B ,|<^ = A^oc. (27) 

As it turns out the parameter A = 1/y/l plays the role of an ultra-violet cut-off A The 
elimination of the matrix element h pq in Eq . fll8|) reminds us to the standard concept of 
renormalization by Wilson, where the high energy modes are integrated out in the path 
integral representation resulting in the effective action for the low energy scales. Indeed, 



3 CANONICAL QED HAMILTONIAN ON THE LIGHT-FRONT 



9 



performing the I (A)-integration to the leading order one gets for the diagonal elements 
hpp,, (for p = p' in Eq.([R|)), 

= E P (A) + SEp(X) — 5E P (A) . (28) 

that defines the connection of energies at different energy scales and coincides with the 
second order of conventional perturbation theory. In the case of QED 3+1 , as the bare 
cut-off A tends to infinity, the sum in Eq. (|28|) diverges, and one has to introduce the 
corresponding counter term. Note that the sum at the upper limit A is regulated by the 
similarity factor in h pq (X), since only the energy differences \E p — E g \ < A are present. The 
ultra-violet renormalization can be attacked with the technique of flow equations order 
by order in a systematic way, and further work is in preparation. 

In the remainder of this paper, the above schematic equations of Hamiltonian flow are 
worked out explicitly for QED 3+1 on the light-cone. 



3 Canonical QED Hamiltonian on the light-front 

Canonical QED3 + i in the front form has been reviewed recently ||. Therefore, only 
the most salient features are recollected in this section, mostly for the purpose to shape 
notation. The Lagrangian density for QED 

C = --F^F^ + + e4-m)ij (29) 

is considered here in the light-cone gauge A + = A + A 3 = 0. Zero modes will be 
disregarded. The constrained degrees of freedom, A~ and ip- (A± = |7°7 ± are projection 
operators, thus ip± = A±ip, and ip = ip + + ipJ) are removed explicitly and produce the 
canonical QED Hamiltonian. It is defined through the independent physical fields A± 
and ip + p2| . To solve the constrained equations for A~ and ip- the auxiliary fields 

A - A - 9 T + 

A+ ~ A+ (^+)2 J > 

$ = ^+ + (mf3 - ia'dxi) ^-^+ , (30) 

are introduced. The fermion current is J M (x) = ^7^. The resulting canonical Hamilto- 
nian H = P + is given as a sum of the free Hamiltonian and the interaction 

H = P + = H + V + W. (31) 

The free Hamiltonian H , the 'kinetic energy', is 

H = \ J dx + d*x ± fi + m2 ^ V±) ^ + H^fA,) . (32) 
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In the 'interaction energy' V + W, the vertex interaction V is the light-cone analogue of 
the minimal coupling interaction in covariant QED and W = W\ + W 2 is the sum of the 
instantaneous-gluon W\ and the instantaneous-fermion interactions W 2 . The latter arise 
from the constraint equations. More explicitly, the interaction is given by 



V 

w 2 



g / dx + d 2 x± J^A 



dx + d 2 x± J* 



9_ 
2 

%- / dx + d 2 x± 



(id* 



r 



id+ 



(fA v * 



(33) 



By definition, the fields \& = \& + + and A^ = [0, A±, A + j are the free solutions which 
in momentum space are parametrized as 



* a (x) 



E 



E 

A 



dp + d 2 p± 
2p+(2vr) 3 
dp + d 2 p± 
2p+(2vr) 3 



(b(p)u a (p, X)e~ ipx + J(p)v a (p, A)e+^) , 
a(p)e>, A)e-*» + at(p)e*(p, \)e"**) . 



(34) 



The Dirac spinors and the polarization vectors are given explicitly in [^]. The single 
particle operators obey the commutation relations 



~a(p),a\p')} = {b(p),b\p')} = {d{p),d){p')} = 5(p+-p + i)5( 2 \p ± -pl)fr 
Inserting the free fields into the Hamiltonian yields for the vertex interaction 

V = g dx + d 2 x± §(x)7 M ^(x)A M (x) 



(35) 



x+=0 

dp\d 2 p 1A [dp£d 2 p ±2 fdptd 2 p ±3 



2p\ 



2pt 



x 

X 



9 £ _ 

sj{2ir) 3 \ lt x 2 ,x 3 J \f2pt 

[{b\p 1 )u a {p l , X 1 )e +i ^ + d{p x )v a {p u X 1 )e^ s 

1% (d\p 2 )vp{p 2 , X 2 )e +i ^ + b(p 2 )up(p 2 , X 2 )e~ ip ^ 
la\ps)el(p3, A 3 )e +i ^ + a(p 3 )e,(p 3 , X 3 )e~^ 



(36) 



The integration over configuration space yields the vertex interaction as a Fock-space 
operator. The integration produces Dirac-delta functions in the single-particle spatial 
momenta, 



dx 



2tt 



d 2 x± 



KT.pI). 
* (2) (E^), 



(37) 
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and often in the sequel these will be used as a 3-dimensional 5-function 

j j j 

The sums over j run over all respective single particles. The Dirac-delta's reflect three- 
momentum conservation, as always in a Hamiltonian approach. W\ and W2 as Fock-space 
operators are obtained correspondingly and found explicitly in M. 

4 Flow equations applied to QED 

In the sequel we consider the canonical Hamiltonian for QED as given in Eq. Q5T| ) and work 
out the details when straightforwardly applying the flow equations as given in Eqs.(|])- 
(lllD . Since we restrict ourselves to solve them up to second order in the coupling constant, 
we can content ourselves to include explicitly only two Fock-space sectors: The sector with 
one electron and one positron (ee) can be identified with the P-space discussed above, 
and the sector with one electron, one positron, and one photon (eej) with the Q-space. 
As a result one aims at the effective Hamiltonian in the ee space. It is helpful to know 
that the effective interaction must turn as the Coulomb potential, to lowest order of 
approximation. 

As a technical trick and for gaining more transparency we omit first the instantaneous 
interactions W. They will be re-installed at the end of the calculation. The physical 
argument is that the instantaneous interaction is already of order g 2 ; the changes due 
to the flow are of higher order in the coupling constant and have to be omitted here 
by consistency. The light-cone Hamiltonian Eq. ( |3T1) is then H = H$ + V and has a very 
simple structure. Since the vertex interaction can not have diagonal matrix elements, both 
diagonal sector Hamiltonians PHP and QHQ are diagonal operators from the outset. The 
case studied here is thus a realization of the paradigmatic case discussed in the second 
part of Sec. [|. To the leading order in the coupling constant the particle number changing 
part PHQ is given by the vertex interaction, i.e. H r (l) = V(l). By reasons to become 
clear soon, we shall put hats on the operators in this section. 

At finite flow parameter / the Hamiltonian is H(l) = Hd{l) + H r (l). The unitary 
transformation of the flow equations diminishes V(l) and generates a new interaction, 
U{1). This interaction is (up to second order) diagonal in particle number and contributes 
to H d {l) = H (l) + U(l). The flow equations Eqs . (|9|) - ([TT|) become then consecutively 

m),H)), (39) 

m),Ho(D\, (40) 
[H (l),V(l)}. (41) 

In this section all of these operators will be evaluated explicitly. The free part is 

Ho{l) = E / dp + d 2 P± E{p; I) (&t(p, \)b(p, A) + d\p, \)d(p, A) + a\p, \)a{p, A)) . (42) 



dU(l) 

dl 
dV{l) 

dl 

rj(l) 
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The single particle energies (E = p ) depend on the 3-momentum p = (p + ,p±) 

m 2 (p; I) + p 2 



E(p-l) 



P 



(43) 



and potentially on the flow parameter through the mass m 2 (p; I) of the particle in question. 
The interaction term is obtained from Eq.(BB) 



V(l) 



(27T) 



X 



[d 3 Pi] f [d 3 p 2 ] f [d 3 p 3 ] {3) 

~T= / ~T= / ~T= 6 (Pi -P2-P2)g(pi,P2,Pa;l) 



2pt 



2p\ 



2p\ 



b\b 2 a 3 (u^ 3 u 2 ) - d{d 2 a 3 (f 2 ^i) + a\d 2 b 3 (v 2 ftu 3 ) + h.c. . 



The integration symbols denote 



[d 3 p] = dp + d 2 p± 



(44) 



(45) 



for example, and the abbreviations u\ = u(p\,\i) and = 7 M e^(p3,A 3 ) are introduced 
for the sake of a compact notation. The effective coupling 'constant' has the initial value 



g(pi,p 2 ,p 3 ;l = 0) =g = e. 



(46) 



with the fine structure constant a = e 2 /4ir ~ 1/137. Correspondingly, the generator of 
the unitary transformation fj = [H , V] becomes 



1 



;2tt) 3 

X 



[d 3 Pi] f[d 3 p 2 ] f[d 3 p 3 



2p+ 



2pl 



2pt 



- h.c. . 



b{b 2 a 3 {uii 3 u 2 ) - d\d 2 a 3 (u 2 j4^i) + a\d 2 b 3 (v 2 ftu 3 )\ - h.c. . (47) 
The structure of fj is very similar to V because H is diagonal, thus 

v(Pi,P2,Ps;l) = g(Pi,P2,P3;l)D(p l ,p 2 ,p 3 ;l) . (48) 
It is convenient to introduce the difference of single particle energies 

D{ Pl ,p 2 , P3 ; I) = E( Pl ; I) - E{p 2 - 1) - E{p 3 - 1) , (49) 
with the £"s being defined in Eq . ([43]) . The derivative dV{l)/dl = [rj(l), Ho(l)] becomes 



dV(Q = _ f [d 3 Pi] r [d 3 p 2 ] f [d 3 p 3 ] §{3) , _ , y(pi,P2,P3-,l) 



dl 



2p{ ■> J2p 



(2*0 



xD(p 1 ,p 2 ,p 3 ; I) b\b 2 a 3 (u^ 3 u 2 ) - d\d 2 a 3 (v 2 fovi) + a\d 2 b 3 {v 2 f x u 3 ) + h.c. 



(50) 



Finally, we calculate the new interactions U(l) which are defined through the derivative 
dll(l)/dl = [rj(l),V(l)]. Their calculation is somewhat cumbersome but straightforward. 
Inserting the six terms from Eq.([0]) and the six terms from Eq. fl4"TD gives 36 terms for U{1) 
which by symmetries reduce to six interactions in the 2-particle sectors. In the present 
work one restricts to calculate U eS , the effective interaction between an electron and a 
positron, as illustrated in Fig. [l|. It has an exchange part and an annihilation part. 
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Figure 2: The graph of the in- 
stantaneous exchange interaction. 
Taken from IL 



Figure 1: The effective interac- 
tion between an electron (e) and 
a positron (e). 



4.1 The exchange part 

The matrix element U ex of the exchange part 

tU»(0 = J[d 3 Pi] J[d 3 p 2 } J[d 3 p[} J[d 3 p' 2 ] 5^(p l +p 2 -p' l -p' 2 ) 

U e x(PuP2;p[,P2,l) b j (pi,\i)b(p' 1 ,\' 1 )d ] (p2,\ 2 )d(p 2 ,\ 2 ) . (51) 
is calculated next to some detail. By direct substitution of Eqs.flHD and (0) one gets 



[d 3 Pl ] f[d 3 p 2 ] f[d s p 3 ] HdVi] f[d 3 p' 2 ] f[d 3 p' 3 



2pl ■> ^2pl ■> J2p? ■> J2p' 2 + ■> J2p 3 



b\bi>a 3 ,al,d 2 d 2 i (' 



dUee;ex(P) _ 1 

dl (2tt)-j ^2pJ 

(S^(p 1 -p' 1 -p 3 )5^(p' 2 -p 2 

(v(Pi,Pi,P3;l)g(p' 2 ,P2,p 3 ;l) + v(P2,P2,P3i l )g(PhPi,P3;l)) 

+ 5 {3) (p' 1 - pi -p 3 )5 {3 \p 2 ~p 2 -p' 3 ) [4^2'a3',4 & l & l'] (fyfov 2) 

{v{Pl,Pl,P3]l)g{P2,P 2 ,P 3 ]l)+V{P2,P 2 ,P 3 ]l)g(p'l,Pl,P3-J))) ■ 



P3) 



(52) 

The commutation relations Eq.(|35"D for the photon induce a three- momentum delta- 
function S^(p 3 — p 3 ). Since 



v(pi,p[,pi -p!;0 

g(p'i,PuP'i - Pi', I) 



-v(p'i,Pi,Pi-Pi,l) 

g(pi,p[:Pi — pi; 



(53) 



all of the integrations in Eq. (|52"D can be performed trivially. The sum over the photon 
helicity is carried out by introducing the polarization tensor p2[ [3| 



(54) 



The null vector 7/ M has the components (r] + , ff±, r]~) = (0, 0, 2) and should not be confused 
with the generator rj. Dropping hence forward the argument I in g or r\ for the reason of 
notational compactness the /-derivative of U^ n becomes 



dl 



'2(2tt) 



0(pi 



Pi ) 



6(p[ 



i + 



Pi 
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{u(p u Xih^uipi, a;)) <v( g ) (t7(^, a' 2 ) 7 ^(p2, a 2 )) 

X 



x (v(p'i,PhPi -pi)g(p2,p' 2 ,P2 -p' 2 ) + g(Pi,Px,p'i -pi)v(p2,P2,P2-p' 2 ))(55) 

The 4-momentum of the photon is denoted by g M , see also Eq.(^) below. The step 
function 6{x) = 1 for x > 0, and zero otherwise. Since #(x) + 9{—x) = 1 the sum of the 
two theta-functions will be replaced by unity in the sequel. This reminds to the calculation 
of the gg-scattering amplitude which was given explicitly in ||. To get the full effective 
exchange interaction one has to include also the instantaneous exchange interaction 



e 



2 



u(pi, Aih^o;, a;)) {v(p' 2 , x 2 )Yv(p 2 , a 2 )) 



u™\i) = - —- v^i^u p^-i/zv-^^ (56) 



To the order considered here it is independent of I, i.e. dU™ st /dl = 0. 
4.2 The annihilation part 

Having been so explicit for the exchange part one can proceed rather quickly for the 
annihilation channel where the calculation proceeds quite correspondingly. One defines 
first the matrix element U an by 

U eS ;an(l) = J [^Pl] J [#P%] J Wl] J W 2 ] ^ (Pi + P2 ~ p\ ~ P 2 ) 

Uan(Pi,P2;p'i,P 2 ,l) b J[ {p 1 ,X 1 )b(p' l ,X' 1 )d\p 2 ,X 2 )d{p' 2 ,X! 2 ) . (57) 
Its /-derivative is defined by the flow equations, i.e. 

d%S»(0 1 ft^i] f[d 3 P2 ] f[d 3 p 3 } f[d 3 p[] r[d 3 p' 2 } r[d%] 



dl J^J .J^J J^J Jty+J J^+J Jy+ 



5 {3) (pi +p 2 - pa) 5 (3 Vi + P 2 - Ps) M4aa, 4' d 2'6i' 



{(v(P3,P2,Pl)g(p 3 ,P2,Pl) +V(P3>P2>Pl) 9(P3,P2,Pl))) ■ (58) 

All integrations can be performed explicitly and one arrives at 

dur = 1 (u( Pl , AO7MP2, a 2 )) djUp) (v( P ' 2 , a^yxk, a;)) 

dl 2(2vr)3 fitifipj P + firtfipY 

x (ri(p 1 +p 2 ,p 2 ,p 1 )g(p[ + p' 2 ,p 2 ,p[) + g(p 1 +p 2 ,p 2 ,p 1 )r](p' 1 +p 2 ,p' 2 ,p' 1 ))(59) 

The 4-momentum of the photon is denoted here by p. The instantaneous interaction in 
the exchange channel 

fjinst^ = e 2 (u(p!, AjjrMgg, A 2 )) (v(p' 2 , X^jM^, X'J) w ^ 



2(2-)3 fi-ffi-+ P +2 
is again independent of / in the lowest non-trivial order of the coupling constant. 
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4.3 Integrating the flow equations for the exchange 

The first order flow equations have been given in Eqs. fl40|) and (|4l|) in operator form. 
After evaluating all matrix elements they reduce simply to two coupled equations: 

J, = -D{pi;p2,P3;l) r]{pi,p2,p 3 ;l) , (61) 

r)(p 1 ,p 2 ,p 3 ;l) = D(p 1 ,p 2 ,p 3 ;l) g(pi,P2,PsJ) • (62) 
Replacing g by the suitably normalized similarity function / according to 

g(pi,P2,Ps; l) = g(o)f(pi,P2,P3, l) = ef(p u p 2 ,p 3 ; , (63) 

see also Eq.(|^), one gets 

f(pi,P2,P3] I) = exp (-1 D 2 (p u p 2 ,p 3 )) (64) 

as the explicit solution. It describes the decay rate of the off-diagonal vertex interaction 
V. However, because of the considerations in Sec. [2] particularly Eq.(P2"D and the consid- 
erations below we want to keep / = f(D;l) as a general function. Correspondingly, we 
rewrite r] as 

1 f d\nf(D;l) \ 
W) = -7^( Ti U(0- ( 65 ) 



D\ dl 

The formal integration of the flow Eqs. ( pTj ) and (^) can be treated more compactly in 
a reasonably short notation. We therefore define the always negative quantities 

D e = pi" - Pi ~ (p[ ~ Pi)~ , 

De = P 2 - p' 2 ~ - (Pa - P2Y ■ ( 66 ) 

They represent the energy differences along the electron and the positron line, respectively, 
and are in simple relationship to both the 4-momentum of the exchanged photon 

% = Pin - Pin ~ = P2n ~ p 2 n ~ Vn~Y > ( 67 ) 

and to the (Feynman-) 4-momentum transfers along the two lines 

Ql = -tfi-Pi) 2 = -q + D e , (68) 
Qt = -( P 2-P2) 2 = -q + D s , (69) 

which need not be equal in a Hamiltonian approach. Since the always positive Feynman- 
momentum transfer Q is a more physical quantity than the energy difference, the D's will 
be substituted in the sequel by the Q's as long as no misunderstanding can arise. In fact, 
we shall use often the mean-square momentum transfer and the mean-square difference 

Q 2 = \{Ql + Ql) = - q -j{D e + D E ), (70) 
SQ 2 = \{Q\ - Ql) = - q Y (D e - D s ) , (71) 
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respectively. The above definitions are also useful to simplify the polarization tensor 
appearing in Eq.(^). Since d liv {q) appears always in combinations with the spinors one 
can make use of the Dirac equation (pi — p[)^u(pi)'y^u(p' 1 ) = and write 

qMPi, AOtMPi, A 'i) = -^fVuVfa, AOt^W, A!). (72) 

One can replace thus in Eq.(|55|) d fiu (q) — > —g^ u + %^(Q/? + ) 2 - With these definitions 
we return now to the problem of integrating up Eq. (p5|) . Substituting Eqs.(^) and (|63|) 
one has for its /-dependent part 

For the formal integration it turns out useful to introduce the abbreviation 

G(D e , D s ) = - Hdl' ^^f^ f{D- e ; /') = ee - , (74) 
Jo dv 

which is not symmetric in the arguments but which satisfies by means of Eq. 

G(D e , D s ) + 0p g , D e ) = & eS + e - e = 1 . (75) 



When /-integrating Eq.(|73|), the /-dependence of D^t) in the denominator can be neglected 
to the order considered here, therefore 

dt'(v(Pi,Pi,Pi -Ptil')9(p2,P2,P2-P2> l ')+v(p2,P2,q\ l')9(Pi,Pi,Q;l')) 

+ TT = ~ e ? T^T + 7^ • ( 76 ) 



o 



The latter combination appears repeatedly, see for example Eq.( |138| ). Putting things 
together the generated interaction Eq.([5f|) becomes 

t]9 en = + e 2 (^ 1 ,A 1 )7^(P / 1 ,A' 1 ))(^' 2 ,A' 2 ) 7 ^(P2,A 2 )) 



x + If 1 > (77) 



and after substituting d^ v (q) 

e 2 (u( Pl , A07^(K, A'J) (v(p' 2 , X' 2 )rv( P 2, A 2 )) 



en 
ex 



2 ( 27 3 SpAJW \T¥?\[zp 



( Q 2 \ (e ee , e e - e \ 
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Adding the instantaneous exchange interaction Eq.(^6|), using Eq.(^), one gets 

e 2 (u( Pl , Ai)y«(pi, AD) {v{p' 2 , \' 2 )Yv(P2, A 2 )) 



u„ 



2(2tt) 



x 



e 



ee , ^ee 



Oee Oee 



for the matrix element of the total exchange interaction. 



(79) 



4.4 Integrating the flow equations for the annihilation 

For the annihilation term we define the (now always positive) energy differences as 

D a = v'x + V'z ~ (Pi + P2Y 
D b = Pi + P 2 ~ (Pi + P2T ■ 



30) 



They are related to the 4-momentum of the photon p M in the t-channel by 

Pm = Pit, + P2 M - V»^Y = Pin + P2n - v»y ■ ( 81 ) 

Rather than the momentum transfer Q the free invariant mass-squares of the initial and 
final states are introduced by 



M 2 a = {p' l +p' 2 f=P + D a 
(pi +P2) 2 = P + D b , 



Ml 



as well as their mean and difference 

1 



M 2 
SM 2 



2 (M 2 a +M b/ 2 



l) = ^(D a + D b ) 



P 



(Ml - Ml) = '-{D a - D b ) , 



(82) 
(83) 



(84) 
(85) 



respectively. The generated annihilation interaction becomes 



-,2 /777 



U 9 



en 
an 



u(px, XiWvfa, A 2 )) (v(p> 2 , A' 2 hWi, AD) 



2(2tt) 



2pTJ2p. 



2p' 2 + J2p' 1 + 



x 



(86) 



The Dirac equation (pi + p 2 ) ^u{pi)^v{j) 2 ) = allows to write 

PfAiipi, Xi)^v(p 2 , A 2 ) = ~^V^(pi, Xi)^v(p 2 , A 2 ) . 



(87) 
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After substituting d au (p) 



-9fiu - r] u r] u M 2 /p +2 one gets for 



en 
an 



,2 Arr 



jjgen 
an 



u( Pl , X 1 )Yv(p 2 , A 2 )) (v(p' 2 , \' 2 h u u( P ' 1} a;)) 



2(2tt) 3 
x [guy + V7„ 



M 2 \ /0 fi 



+ 



e 



ba 



Ml 



Adding the instantaneous term yields the effective interaction in the annihilation channel 



(u{ Pl , A 1 ) 7 ^(p 2 , a 2 )) (v{ P ' 2 , KhMti, a;)) 



2(2tt) 



9nv 



e 



ab 



M 2 



+ 



e 6a \ 



e 



M 2 , 



ab 



2p' 2 + J2p' + 



e 



ba 



1 

SM 2 



Ml 



Ml 



P 



39) 



all in perfect analogy to the exchange term. 
4.5 The effective ee-interaction 

Thus far, we have been studying the structure of the Hamiltonian proper H = P + = P /2. 
In dealing with its spectra it is advantageous to study the spectrum of the 'light-cone 
Hamiltonian' H^c = P n P^ = P + P~ — P 2 . It is a Lorentz scalar || with the eigenvalues 
having the dimension of an invariant mass-squared. Combining the effective interaction 
of the exchange and annihilation channels one introduces therefore 



U, 



cff 



p + \u ex + u a 



(90) 



see Eq.(|79D and Eq.fl89|). One gets 



U, 



a 



cff 



4tt 2 
a 
'i^ 2 



< YY >ex 



< YY >an 



9nv 



We 

mi 



+ 



" Ql 

Oba 



Ml 





6Q 2 




e e - e 


V 2 


W 


Ql 




6M 2 


(Oab 


®ba 




p+ 2 


w 


Ml 



(91) 



The symbols < 7^7^ > are introduced conveniently as 

(u( Pl , Ai)7"w(pi, A'J) (v(p' 2 , \' 2 )rv( P 2, A 2 )) 



< YY > ex -- 



< YY >an~ 



{YJ V ) ( 



>xx'(l — x)(l — x') 
(u( Pl , A 1 ) 7 ^ (P2, A 2 )) (v(p' 2 , X' 2 )Yu( P [, AD) 
xx'(l — x)(l — x') 



xx'{\ — x)(l — x') 

V Jan 

xx'{l — x)(l — X 1 ) 



(92) 



where for example x = pf / P + is the longitudinal momentum fraction 
The effective interaction U e g is the kernel of the integral equation 



M 2 (x, k±; Ai, A 2 



m 2 + kf 



(x, k±; Ai, A 2 



+ 



x(l — x) 

/ dx'd 2 k'j_ (x, k±; Ai, X 2 \U e s\x', k\; X[, X' 2 ) (x' , k\] X\, X' 2 \ 



(93) 
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In the equation appear only intrinsic transversal momenta k± and longitudinal momentum 
fractions x = pf/P + , defined by i.e. Pi = (xP + ,xP± + k±,pT). Its spectrum is thus 
manifestly independent of the kinematical state of the bounded system, particularly of 
P + and P±, which reflects the boost invariances peculiar to the front form ||. The 
integral equation replaces in some way Eq. ([!]). The first term on the r.h.s is the free part 
of the Hamiltonian in analogy to a 'kinetic energy', and the second term is an 'interaction 
energy' which is the relativistically correct interaction, correct up to the second order in 
the coupling constant. 

The integration domain D is restricted by the covariant cut-off condition of Brodsky 



and Lepage [22 



1 \Z \ A im 2 



< A + Am 2 , (94) 



x(l — x) 

which allows for states having a kinetic energy below the bare cut-off A. 



4.6 Dependence on the cut-off function 

Before discussing the dependence of the effective potential in Eq. fl9~I|) on the cut-off func- 
tion / one has to determine the dependence of on its two energy arguments D, as given 
in Eq. (|74|) . It is natural to assume that the similarity function f(D; I) is a homogeneous 
function of its arguments 

f(D;l) = f(D K l) (95) 

with some exponent k. As examples we consider three types of similarity functions and 
the corresponding functions: 
the exponential cut-off 

/(£>;0=exp(-|£>|0, « = 1, G(D e , D s ) = — ^— , (96) 



D K + D, 



the Gaussian cut-off 



f(D;l) = ex V (-D 2 l), K = 2, Q(D e , D s ) = — ^— , (97) 



D 2 + D 2 



and the sharp cut-off 



/(£>; /) = 9(1/1 - \D\") = 9(1 - \D\ K l), S(D e , D s ) = 9(\D e \ - \D g \) = 9(1 - ^). (98) 

In the last case k is an arbitrary positive number. In the first and in the last case we have 
used that D e and D s have the same sign. The second case corresponds to the solutions 
of the original flow equations, see Eqs.([L8|) and fl64|) . 

Assuming this homogeneity of /, is quite generally a function of the ratio of its two 
arguments since 
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Instead we may express it as a function of 

£ = 
e(D e ;De) = 



De-D-e 
D P + D- P 



5Q 2 



Q 2 



2 (* + *(£))• 



(100) 
(101) 



Due to Eq.(^D is an odd function, #(-£) = -0(f)- 

Then the first expression in square brackets in ([Jl]) which contains the singular part 
of the effective interaction reads 



rjex 




— ec — cc i & / — cc — cc 

7^2" TyT ^p^ v ~~T2 ~7y2 7^7 



Q 2 



^g 2 




V 2 








q +2 J 





q 2 i-e 

For the three similarity functions mentioned above one obtains 

e - zm 



exponential = £, 

Gaussian = 3^2, 

sharp tf(f) = sign(f), 



e 2 -^(o . -e 2 
e 2 -^(o -lei 



(102) 

(103) 
(104) 
(105) 



We observe that the effective interaction depends explicitly on the similarity function. The 
requirement of block diagonalization of the Hamiltonian determines the generator only 
up to a unitary transformation of the blocks. This explains why the effective interaction 
may depend on the similarity function. 

We will discuss the dependence on this function further in the next section, but mention 
that for the elimination of the electron-phonon interaction in solid-state physics which 
yields the effective attractive interaction between electron pairs responsible for super- 
conductivity, one may also choose different similarity functions, see |12| and | 



21 . For 



realistic spectra, Mielke |2T| has found that the critical temperature calculated from the 
Gaussian similarity function and that suggested by Glazek and Wilson Q differ by only 
2%, the difference to those calculated by the conventional Eliashberg theory was only 5%. 

Since the kernel of the integral equation is manifestly frame-independent, one can 
evaluate it in the particular frame P + = 2m and P± = 0. For the further discussions we 
choose to express the momenta as 



Pi 

pi 



m + 
m — 



and similarly for p[ and p' 2 . With g H 
energy differences become 



P\\, Pi± 

P\\, P2± 

Q\\ = P'\\ 



Px 
-P± 



(106) 

p 2 = p]_ + p| and q 2 = qj_ + q 2 the 



m 2 + p' ± 



m 



PL 



m + 



P|| 



m + 



ql 

q+ 



2 2 /2 

q p p 

gy m — p\\ m — pj, ' 



(107) 
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m 2 +pl m 2 +p'l ql q 2 p' 2 p 2 

V s = = j H . 

m — p|| m — pj| q + q\\ m — p|| m — py 
The other quantities become correspondingly 

D e -D g = 2ml / /2 - f I , (109) 

= (110) 
Q 2 = 5 2 + si f — r^— ^ — ^"72) ■ (m) 

ym 2 -p|| 2 m 2 -p{ J 

Note that — m < pii < m and that pn may not be interpreted as the z-component of a 
single-particle momentum. 



5 Discussion and interpretation 



The integral equation (|93|) seems to have two kinds of singularities: The 'Coulomb singu- 
larities' 1/Qg or l/Q? and the 'collinear singularity' Either denominator can become 
zero in the integral equation. The Coulomb singularity is square-integrable and welcome 
since it provides the binding. The collinear singularity is disastrous. If the coefficient of 
l/q?< is finite at q\\ — > the integral equation is not solvable. The singularity structure of 
the final result must therefore be discussed carefully. 

First we observe, that in the case of the exponential similarity function the effective 
interaction between an electron and a positron becomes 

4?r 2 Jx(1-x)x'(1-x') V Q 2 M 2 J 



The collinear singularity is wiped out since the coefficient related to the similarity function 
vanishes, see Eq.( |103| ). This astoundingly compact formula exactly agrees with the Tamm- 
Dancoff approach ||10|| . The explicit x-dependence in the denominator of Eq.(|9~2"l) looks like 
the only remnant of the light-cone formulation; all other quantities are Lorentz-contracted 
scalars. 

With the method of Hamiltonian flow one can calculate also the scattering amplitude 
for ee-scattering 

rp a 1 ( < j% > ex < c^C^ > an \ 

SCattermS " WJxd-xWl-M Q 2 M 2 J 1 {U6} 



sec App. ^ The expression agrees identically with the Feynman amplitude. In the 
scattering process the free 4-momentum P M is conserved, and thus the momentum transfer 
along the electron and the positron line are the same, i.e. Q\ — Q\ — Q 2 ■ Therefore 
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SQ 2 = and thus £ = 0. According to Eq. ( |102|) , the coefficient of the collinear singularity 
vanishes identically. 

The scattering amplitude can be investigated also in light-cone perturbation theory 



. Its exchange part is presented in all details in Sec. 3.4 of 0. As can be pursued 
there, the collinear singularity appears in both the instantaneous interaction and the 
second order amplitude VGV of the vertex interaction, however such that the two contri- 
butions cancel each other exactly. Only the integrable Coulomb singularity remains. One 
concludes that Hamiltonian flow equations are governed de facto by the same mechanism 
of cancelation but that the disappearance of the collinear singularity is somewhat more 
subtle. To understand that better two approximations are discussed in the sequel. 

First, let us expand the kernel up to terms linear in p/m. To this order holds 

u(p', aOtMp, a) = { 2m ^ ^ZX] I' ( 114 ) 

The bracket symbols of Eq. ( j92|) become then 

< Yl u > =< Yl u > Vtfu = 16m 2 <$<$.. (115) 

According to Eqs. ([L07|) and (|108|) one has to this order 

D e = D- e = -q 2 /q\ V (116) 

Consequently, 5Q 2 = and that alone is sufficient to wipe out all dependence on the 
similarity function. The exchange part becomes 

U eS = -~(2m) 2 ^. (117) 

The factor (2m) 2 is a light-cone peculiarity, see Sec. 4.9 of ||, and the factor in front 
of it is precisely the Fourier transform of the familiar Coulomb potential in three space- 
dimensions. Hence the effective electron-positron interaction Eq.(|9"lf) is bound to produce 
the Bohr spectrum. 

Next, expand the kernel up to second order in p 2 /m 2 . This gives the familiar Breit- 
Fermi spin-spin and tensor interactions ||, that insures the correct spin-splittings for 
the positronium ground state and restores the rotational invariance J7|. Brisudova et al. 
H state also that the 'rj^rj^ term may influence the spin-orbit coupling in second-order 
bound-state perturbation theory. The correct singlet-triplet splitting is observed also in 
the numerical solutions of the integral equation [TU, |T7| with the effective interaction ( |112j ) 
as well as with flow equations (see further). 

Let us now consider the case of a general similarity function. In order to do this 
we have to discuss the expressions of the previous subsection. We realize, that as the 
momentum transfer q = p' — p tends to zero D e — D s tends linearly to zero and Q 2 
vanishes quadratically in q. This approach to zero is anisotropic. Only at p = itself 
Q 2 approaches isotropically and the linear contribution in D e — D s disappears. One 
can also show, that Q 2 can only vanish, if q — 0. It is not sufficient, that q + vanishes. 
Consequently £ is always finite (which also follows from the observation, that Q 2 e and Q 2 
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are positive and thus — 1 < £ < +1). However as q vanishes £ will in general depend on the 
direction from which q approaches zero. Moreover we realize, that for sufficiently smooth 
similarity functions / the function will be analytic as in the cases of the exponential 
and Gaussian cut-off. Then the pre- factor (£ 2 — /(l — £ 2 ) in ( p. 02] ) contains a factor 

£ 2 , which itself contains a factor q +2 , which cancels against the denominator q +2 of the 
rj^rjy term. Thus the interaction becomes only singular if q approaches where it diverges 
like 1/q 2 . This is however not true for the sharp cut-off, where only one factor q + can be 
cancelled and one q + remains in the denominator of the t]^r] u term. Thus for a smooth 
cut-off one gets rid of the collinear singularity. 

We realize further that £ is of order p/m. Since however the Bohr momentum is of 
order ma the contribution due to the second term in ( |102[ ) is smaller by a factor a 2 in 
comparison to the leading term g^/Q 2 . 

Thus we realize that in this order we have in addition to the leading term a contribution 

- £i < tv >« (i - m)e (ff - ^) (us) 

to the interaction. It obviously depends on the similarity function via the derivative $'(0). 



We will shortly discuss the leading contribution and use Eq. (|115|) . Collecting terms one 
obtains 

(l _ m) ]^<Wl_f? (119) 

as the contribution to the effective interaction. This interaction is spin-independent but 
anisotropic. It depends on the similarity function and is of order a 2 in comparison to 



that of the leading Coulomb interaction. We emphasize as already pointed out in JL3 
that in this order in a 2 also one- and two-loop terms contribute. Since they will also in 
general depend on the cut-off chosen, there should be a cancelation between the term 
found here and loop terms. It was argued there, that in order a 2 the one- and two-loop 
terms only contribute spin-independent interactions and interactions of one spin (spin- 
orbit coupling), but not interactions between both spins (spin-spin and tensor interaction), 
which is in agreement with our finding that the interaction ( |119| ) does not depend on the 
spin. 



6 Numerical solution for positronium spectrum 

We solve the integral equation Eq.(^), with interaction kernel given in Eq. (|9~TD , for 
positronium mass spectrum numerically. Effective interaction with different choice of 
cut-offs is sumarized in Appendix (|FJ). 

6.1 Formulation of the problem 

In polar coordinates the light-front variables are (k±;x) = (k±,(p;x); therefore the ma- 
trix elements of the effective interaction Eq.(|9T|) depend on the angles (p and (p f , i.e. 
(x, k±, if] Ai, A2|Vrff|x / , k' ± , <f'; X[, X' 2 ). In order to introduce the spectroscopic notation for 
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positronium mass spectrum we integrate out the angular degree of freedom, <p, introduc- 
ing a discrete quantum number J z = n, n e Z (actually for the annihilation channel only 
| J z \ < 1 is possible), 



(x, k±; J z , Ai, A 2 |V r efr |a; / , k' ± ; J' z , A' 1; A' 2 ) 

= -I 

2tt Jo 



' 1 ^ d<fe- iLz(p d(p'e iL *v'(x,k ± ,<p; X 1 ,X 2 \V eS \x',k' ± ,(f'] X[,X' 2 ) 



(120) 



where L z — J z — S z ; S z — ^-+ 2 anc ^ the states can be classified (strictly speaking only for 
rotationally invariant systems) according to their quantum numbers of total angular mo- 
mentum J, orbit angular momentum L, and total spin S. Definition of angular momentum 
operators in light-front dynamics is problematic because they include interactions. 

The matrix elements of the effective interaction before integrating over the angles, 
(x,k±,ip; Xi, X 2 \V e g\x' ,k' ± ,<p'; X[, X 2 ) , and after the integration inroducing the total mo- 
mentum, J z , (x, kj_; J z , Ai, A2I V^ffl^;', k' ± ; J' z , X[, X' 2 ) for different cut-off functions are given 
in the exchange and annihilation channels in Appendices (|C]) and (|B|), respectively. 

Now we proceed to solve for the positronium spectrum in all sectors of J z . For this 
purpose we formulate the light-front integral equation Eq. (|9~3"D in the form where the 
integral kernel is given by the effective interaction for the total momentum J z , Eq. ( |120| ). 
We introduce instead of Jacobi momentum (x, k±) the three momentum in the center of 
mass frame p = (p z , k±) = (fi cos 9, \i sin 6 cos (p, fi sin 6 sin <p) as follows 

(121) 

where the Jacobian of this transformation dx/dp z is 

_ 1 m 2 + fi 2 sin 2 6 

~ 2 (m 2 + ^ 2 )3/2 ' [LZZ) 

One obtaines then the integral equation 

(M 2 - 4(m 2 + /i 2 ))^„(/x, cos^; J z , X h A 2 ) 

1 n/ jj' 2 m 2 + /i' 2 (l - cos 2 6>') 



+ y [ du! [ dcosO' 

T , w M JD J-l 



2 (m 2 + fi' 2 ) 3 / 2 

x (ft, cos6*; J z , Ai, A2 1 Vefr I A*', cos^'; J' z , X[, X 2 )ip n ({i', cos9'; J' z , X[, X 2 ) = 0. (123) 



The integration domain D, defined in Eq. (|94]) , is given now by fi e [0;f]. Neither L z nor 
S z are good quantum numbers; therefore we set L z = J z — S z . 

The integral equation Eq. ( |123|) is used to calculate positronium mass spectrum nu- 
merically. Note, that if one succeeds to integrate out the angular degrees of freedom 
for the effective interaction Eq. ( |120| ) analytically, one has 2-dimensional integration in 
Eq. ( |123| ) instead of 3-dimensional one in the original integral equation (|93|) to perform 
numerically. 
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n 


Term 


Betpt 


B E 


Bq 


B G 


B s 


1 




1 118125 


1 049550 


1 101027 

X • X \J X UiJ I 


1 026170 


920921 


2 




0.998125 


1.001010 


1.049700 


0.981969 


0.885347 


3 


2% 


0.268633 


0.260237 


0.266490 


0.260642 


0.242607 


4 


2 3 ^ 


0.253633 


0.253804 


0.259506 


0.254765 


0.234312 


5 


2 1 P 1 


0.253633 


0.257969 


0.263056 


0.257664 


0.237611 


6 


2 3 P 


0.261133 


0.267070 


0.273826 


0.266563 


0.243075 


7 


2 3 Px 


0.255508 


0.259667 


0.265412 


0.260127 


0.238135 


8 


2 3 P 2 


0.251008 


0.255258 


0.260345 


0.255498 


0.236383 



Table 1: Binding coefficients, B n = 4(2 — M n )/a 2 (a = 0.3), for the lowest modes of the 
positronium spectrum at J z = for the equal time perturbation theory up to order a 4 
(Betpt |10| ) compared to our calculations with exponential (Be), Gaussian (Bq) and 
sharp (B s ) cut-offs. B G is obtained using only g^ v part of interaction; for B G 't]^t] u ' term 
is included. Exchange channel is considered. 



We use the numerical code |ilf , worked out by Uwe Trittmann for the similar problem 
This code includes for the numerical integration the Gauss-Legendre algorithm 
(Gaussian quadratures). To improve the numerical convergence the technique of Coulomb 
counterterms is included. The problem has been solved for all components of the total 
angular momentum, J z . 

Positronium spectrum is mainly defined by the Coulomb singularity 

g-0, (124) 

which is an integrable one analytically and also, by use of technique of Coulomb countert- 
erms, numerically. In this region the effective interaction Eq.(|9~T|) has leading Coulomb 
behavior Eq.( |112j ), independent on the cut-off function. We use therefore in numerical 
procedure standard Coulomb counterterms, introduced for the Coulomb problem Eq. ( |112| ) 
0,0, for all cut-offs. 

Also we expect therefore the same pattern of levels for different cut-offs, that is proved 
numerically further. 

Another important limiting case to study effective interaction Eq . (|9l|) , namely its 
exchange part Eq. ( |102j ), is the collinear limit 

g+^0, (125) 

that is special for light-front calculations. Exchange part of the effective interaction is 
given by Eq.( |118| ), which is finite in this limit. This is true for the regular cut-off func- 



tions, as in the case of exponential and gaussian cut-offs, where the derivative d'd(0)/d£ i 
is well defined. For sharp cut-off this condition is not fulfilled, and the effective interac- 
tion contains the l/q + type of singularity (see Appendix (||). We do not associate any 
physics with this singularity, and consider it as a consequence of artificial choice of cut-off, 
which corresponds to singular generator of unitary transformation Eq.(^). In numerical 
calculations we omit 'rj^rjl term in exchange channel for sharp cut-off. 
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n 


Term 


B E 


Bl 


B G 


B s 


1 


1% 


1.049550 


1.101270 


1.026170 


0.920921 


2 


l 3 ^ 


0.936800 


0.978018 


0.921847 


0.834004 


3 


2% 


0.260237 


0.266490 


0.260642 


0.242624 


4 


2 3 S 1 


0.255292 


0.260383 


0.255615 


0.234338 


5 


2 1 P 1 


0.257969 


0.263056 


0.257664 


0.236383 


6 


2 3 P 


0.267090 


0.273847 


0.266626 


0.243075 


7 


2 3 Px 


0.259667 


0.265412 


0.260127 


0.237611 


8 


2 3 P 2 


0.245615 


0.250821 


0.247091 


0.230901 



Table 2: Binding coefficients, B n = 4(2 — M n )/a 2 (a = 0.3), for the lowest modes of the 
positronium spectrum at J z = for our calculations with exponential (Be), Gaussian 
(Bq) and sharp (Bs) cut-offs. B G includes 'rj^rj^ term in exchange channel; Bg does not. 
Exchange and annihilation channels are considered. 



n 


Term 


1 SB E 




1 SB G 




I SB S 




2 


l 3 ^ 


6.30 10- 


4 


1.76 10" 


-3 


1.18 10 


-3 


4 


2 3 S' 1 


8.40 10- 


5 


1.77 10" 


-4 


9.0 10- 


5 


5 


2 1 P 1 


-1.30 10" 


-5 


-7.47 10 


-4 


-9.1 10" 


-5 


7 


2 3 P X 


-4.08 10" 


-4 


-4.08 10 


-4 


1.4 10- 


-4 


8 


2 3 P 2 


5 10- 6 




-7.7 10- 


5 


4.15 10 


-4 



Table 3: Difference in the corresponding energy levels between J z = and J z = l states 
for exponential (5 Be), Gaussian (SB a) and sharp (SBs) cut-offs. Exchange channel is 
considered. 



n 


Term 


j 5B E \ 


SB G 




I 5B S 




2 


l 3 ^ 


-1.411 lO" 3 


-7.86 10- 


-4 


-1.65 10- 


-3 


4 


2^ 


-4.1 10~ 5 


-4.0 10- 


5 


-1.15 10- 


^4 


5 


2 1 P 1 


-6.4 10~ 5 


-6.52 10- 


-4 


-4.60 10- 


-4 


7 


2 3 P X 


-4.69 10" 4 


-4.74 10- 


-4 


-1.40 10- 


-4 


8 


2 3 P 2 


-1.96 10~ 4 


-1.36 10" 


-4 


-2.44 10- 


-4 



Table 4: Difference in the corresponding energy levels between J z = and J z — 1 states for 
exponential (8 Be), Gaussian (8Bq) and sharp (SBs) cut-offs. Exchange and annihilation 
channels are considered. 
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6.2 Discussion of numerical results 

We place the results of calculations for three different cut-offs, performed in exchange and 
including both exchange and annihilation channels, in Tables ([[]) and (0), respectively. 
The corresponding set of figures is presented in Fig.(|3|) and Fig.(^). We get the ionization 
threshold at M 2 ~ 4m 2 , the Bohr spectrum, and the fine structure. Including annihilation 
part increases the splittings twice as large for the lowest multiplets. 

We argued that the region of Coulomb singularity, and hence 'g' part of effective 
interaction, determines mainly the positronium spectrum. However, including 'rj^l part 
for gaussian cut-off shifts spectrum as a whole down to about 5 — 7%, since this part 
is diagonal in spin space (Appendix (0)), and improves the data to be near the result 
obtained in covariant equal time calculations (Table (|l|)). For the sharp cut-off the lowest 
multiplet is placed higher than the one in case of exponential and gaussian cut-offs. The 
reason is in disregarding the infrared divergent 'rj^rjl part. Presumably, it is necessary to 
take into account 'rj^rjl term in exchange channel also for sharp cut-off after the proper 
regularization of infrared longitudinal divergences is done. 

As one can see from presents figures, certain mass eigenvalues at J z = are degenerate 
with certain eigenvalues at other J z to a very high degree of numerical precision. As an 
example, consider the second lowest eigenvalue for J z = 0. It is degenerate with the lowest 
eigenvalue for J z = ±1, and can thus be classified as a member of the triplet with J = 1. 
Correspondingly, the lowest eigenvalue for J z = having no companion can be classified 
as the singlet state with J = 0. Quite in general one can interpret degenerate multiplets 
as members of a state with total angular momentum J = 2J z ^ max + 1. One can get the 
quantum number of total angular momentum J from the number of degenerate states for 
a fixed eigenvalue M 2 . One can make contact with the conventional classification scheme 
2S+1 L J f, as indicated in Tables ©-(0). 

Such pattern of spectrum is driven by rotational invariance. To trace rotational sym- 
metry we calculate the difference of energy levels between J z = and J z —\ states for the 
lowest multiplets. The data are given for exchange and including annihilation channnel 
in Tables (||) and @, respectively. Including annihilation channel improves the extent of 
degeneracy (see Table (f|) and Figure^)). 

Concerning the spin-splittings the best agreement with covariant calculations is ob- 
tained for gaussian cut-off, the worst results are for sharp cut-off. Rotational invariance is 
traced on the level of spectrum by studing the degree of degeneracy of corresponding states 
with the same total momentum but different projection J z in the multiplet. Again, better 
results are obtained for exponential and gaussian cut-off functions than for sharp cut-off. 
This suggests, that smooth cut-off functions are preferable to perform calculations. 

Generally, the impact of the different choice of cut-off functions on the spectrum is 
small. 

In this work we solve the bound state integral equation for the one fixed integration in- 
terval. Integration domain introduces the ultraviolet cut-off dependence of invariant mass 
squared M 2 (A), that reflects renormalization group properties of the effective coupling 
constant. We leave this question for the future study. 
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Figure 3: The invariant mass-squared spectrum Mf for positronium versus the projection 
of the total spin, J z , excluding annihilation with exponential, Gaussian and sharp cut-offs. 
The number of integration points is N ± — 7V 2 = 21. 
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Figure 4: The invariant mass-squared spectrum Mf for positronium versus the projection 
of the total spin, J z , including annihilation with exponential, Gaussian and sharp cut-offs. 
The number of integration points is N± = N 2 — 21. 



7 Summary and conclusions 

We have applied the method of Hamiltonian flow to the canonical Hamiltonian of quantum 
electrodynamics in the front form and derived an effective interaction for an electron and a 
positron, which acts only in the ee-sector of Fock space. To lowest order of approximation 
is the familiar Coulomb interaction. In this first of a series of papers we have restricted 
ourselves to include terms up to the second order in the coupling constant e. By reasons of 
simplicity (almost) all aspects of renormalization theory within a Hamiltonian approach 
have been disregarded in this first assault. Depending on the particular choice of the 
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similarity function one gets perfect agreement with other approaches particularly with 
the method of iterated resolvents. Special emphasis is put on the impact of the collinear 
singularity. Depending on the similarity function, the collinear singularity is either absent 
from the outset, or it is shown explicitly that it has no impact on the solubility of the 
final integral equation. 

The numerical solution of positronium bound state problem, with the effective electron 
positron interaction obtained by the flow equations, is presented. No approximations along 
numerical procedure are done. 

One concludes that the method of Hamiltonian flow equations looks like an excellent 
tool to progress further, particularly to attack the severe problems of renormalization 
theory within a non-perturbative Hamiltonian approach to field theory and thus to apply 
it eventually to non-abelian gauge field theory and quantum chromodynamics. 
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A Flow equations and Feynman amplitudes 

There must be a connection between the Feynman diagrammatic technique which is based 
on the action and the flow equations which are based on the Hamiltonian. It will be studied 
here up to second order in the coupling constant 

The difference resides in the /-ordering of the generator r)(l) in the unitary operator 

U(l) = Ttexp ^ r](l')dl^j . (126) 

The transformed Hamiltonian reads 

H{1) = U(l)H{0)U + (l) . (127) 
With the full unitary transformation 

U(oo) = e~ s (128) 
one has thus up to the second order in the coupling constant 

H(oo) = e- s H(0)e s = H(0) + [H(0), S] + [[if (0), S], S] + ... , (129) 

and 

S = + S {2) + ... = -/ dl rj(l) - - / dl / dl' [t/(0, r]{l')\ + ... . (130) 

jo 2 jo jo 

The series for the effective Hamiltonian reads 

fi(oo) = H (0) + V(0) + [H (0),S^] + [V(0),S^} 

+ \[[H o (0),S^%S^] + [H o (0),S^] + .... (131) 

The first order term vanishes 

V(P) + [H Q (0),SW] = 0, (132) 

which gives rise to 

if (oo) = H o (0) + I[K(0), S™] + [ff o (0), S (2) ] . (133) 
Denoting by the order with respect to e and omitting the instantaneous terms gives 

= ^/l^i/^l/l^l^)(p 1 -p 2 -P3) 1 



X 



b\b 2 a 3 (uifou 2 ) - d{d 2 a 3 (u 2 ^i) + a[d 2 b 3 (v 2 <ftu 3 ) - h.c. . (134) 
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and 



E 



Ai,A 2 ,Ai,A£ 



[d*pi] f[d s p 2 ] r[d%] r[d%] 



2p'{ 



1+ 



2p\ 



2 + 



x 8®{ Pl + P2 - pi - P ' 2 ) (-— J^) (ep X) a) - ep 2) D] 



)) 



x 



(?t(pi, Ai)r«(pi, a;)) w 2 , a 2 ) 7 ^(p2, a 2 )) 



x :tf(p 1 ,\ 1 )b(p' 1 ,\' 1 )J(p 2 ,\ 2 )d(p 2 ,\ 2 ) : 



(135) 



The .D are of course the energy differences defined in Eq. (|66|) , particularly D\ = p[ — 
Pi — (p[ — pi)~ and D 2 = p' 2 ~ — p 2 — (j> 2 — P 2 )~, and q denotes the 4-momentum of the 
exchanged gluon. The matrix elements in the ee-sector become correspondingly 



< ee\[V{0),S (1) }\e'e> > 



{u{ Pl , AxhMPi, A;)) (v(p' 2 , \ 2 )Yv{p 2 , A 2 )) 



x 



2(2vr) 3 



1 



1 



(136) 



and 



<ee\[H (0),S {2) ]\e'e' > 



,2 (tt; 



u( Pl , x 1 )ru(p[, a;)) (v( p ' 2 , x' 2 )rv( P2 , a 2 )) 



x 



2(2vr) 3 
dnv{q) 



'2ptpp+ 

ID 1 -D 2 



2 D X D 2 



2p 2 + \f2pl 
(e(D 1 ,D 2 )-e(D 2 ,D 1 )) . (137) 



Inserting the latter two equation into Eq. (|133|) gives the generated interaction in Eq.(|7 
Note that the 0-factor in Eq.(77) satisfies 



ei 2 e 21 



Do 



-{- 

2 \D 1 



1 

D~ 2 



1 
~D~i 



1 

~~D 2 



(e 



12 



e 



21, 



(138) 



The first term in Eq.([77|) corresponds therefore to the result of perturbation theory and 
comes from Eq.(136). The second term originates from the Z-ordering ( 137 ) and vanishes 
on the mass shell. 
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B Defining different cut-offs 

In this appendix we summarize the results for the effective electron-positron interaction, 
generated by the flow equations with different similarity functions. In the practical work, 
three different similarity function will be studied explicitly: 

(1) the exponential cut-off, (2) the gaussian cut-off, and (3) the sharp cut-off. 
(1) Exponential cut-off 



f(D;l) 



exp 



\D\l) 



D f 



a .( 7 Y) 9 ^ 



471-2 



q+D 
1 

Q 2 



(7Y) £ 



a 

47T 2 



9iw_ 
p+D 
1 

M 2 ' 



(139) 



where D = l/2(D e + D s ) and D = l/2(D a + D^). The first choice of similarity function 
gives exactly the result of perturbation theory. 
(2) Gaussian cut-off 



f(D;l) 
&(D e ,D s ) 



exp 



-D I 



Dl 



D 2 + Dl 



cff 



a 

4^ 
a 

a 

a 
Arc 2 



g^u D e + D- e 



VfiVu (-D e - D e 



q+ Dl + Dl 2q+ 2 D 2 e 

g^v Da + Db _ VyVv (Da - Db) 2 
p+Dl + Dl 2p+ 2 Dl + Dl 



9^ VnVv $Q 4 



Q 2 q+ 2 Q\ 

9ixu _ V7^ SM 
M 2 



Q 4 



v 



M 4 



Q 4 + 5Q± 
M 4 
M 4 + 5M 4 



(140) 



where we understand under Q 4 = (Q 2 ) 2 and 5Q 4 = (5Q 2 ) 2 with Q 2 and 5Q 2 defined in 
Eq. flTTD. 

(3) Sharp cut-off 



f(D;l) 
®(D e ,D s ) 



Vr 



0(1- 

9{\D 
a 



\D\l) 

- \D S \) ; m = sign(e) 



cff 



471-2 
V7" 



2q 



2 (D e -D 



9iw_ 
q + 



\D e \-\De\) , 9(\D s \-\D e \Y 



\De\ ~ \D e \) 



D r 



\D, 



\De\Y 



D P 



D, 



a 
An 2 



(tV>c 



P + { 



\D b \) 



\D h 



\Da\Y 



D h 
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2p+ 2 



(A, - D 



a 
Air 2 

a 
Air 2 



(7Y) £ 




I A 



\Da\Y 



D h 



9nu V»Vv \5Q 

Q2 + „+2 r>2 



M 2 



Q 2 



p 



M 2 



Q 2 + \5Q 2 \ 
M 2 
M 2 + \SM 2 \ 



The motivation to choose these cut-off functions is the following. Using exponential cut- 
off in flow equations one generates the same interaction as obtained also in Tamm-Dancoff 
approach, where numerical calculations of positronium spectrum are performed ||10|| , and 
we use this numerical code here. Note also, that for this cut-off the effective interaction 
looks very much as in covariant calculations: it contains only ' part, and 'rj^rjy part 
is identically zero, so that there is no collinear problem. Gaussian cut-off corresponds to 
the original choice of generator Eqs.(|5|,f4TD by Wegner as commutator of diagonal, particle 
number conserving, and off-diagonal, particle number changing, parts of Hamiltonian. 
Sharp cut-off is used often in the alternative similarity scheme to perform calculations M . 
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C The matrix elements in the exchange channel 

In this Appendix we follow the scheme of the work [TIJ to calculate the matrix elements 
of the effective interaction in the exchange channel. 1 Here, we list the general, angle- 
dependent matrix elements defining the effective interaction in the exchange channel and 
the corresponding matrix elements of the effective interaction for arbitrary J z , after in- 
tegrating out the angles. Exchange part of the effective interaction for three different 
cut-offs Eqs. ( 139 -141) can be written 



Of. 



where explicitly one has 
(1) Exponential cut-off 



(2) Gaussian cut-off 



V ea = - — (Yl u )B, U} (142) 



B„ = %, (143) 



Q' 2 



= ^ Re [WTWi) ~ ^ 5QHm (qw) ' (144) 

where Re and Im are real and imaginary parts, respectively, and i 2 = — 1. 
(3) Sharp cut-off 

(6i-m_ 6(5Q 2 ) \ 



where q = p[ — px is the momentum transfer; and (7^7^) for the exchange channel is given 
in Eq. (p2|). We omit index 'ex' everywhere. 

It is convenient to extract the angular dependence in the functions 

Ql = a\ — b cost 
Q\ = &2 — boost 
t = <p-<p', (146) 

where we define 

k± = k^(cosip,sm{p) (147) 
in polar system; here the terms are given 



ry>' nn I ry> /y' 

JU ^ r) Jj ^ / <2 \ ,jb ,Xj 



l\1 



ax = — k ± H — -kf + m 



x x' XX 1 



1 Some of these calculations can be found in mA 
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= k l + k 'l + ( x - x ')(kl(--)-k'l(-^)]+m 2ix X>) 

\ x x J 

rf ^ 1 rf ( rf rf^ \ 2 

Jb ty ± Jb fry ry IX tX/J 

a 2 = fcl + -k_r + m 2 



l\2 



xx' 



1-x x 1-x' (1 -a;)(l -a;') 

fc 2 + + , x _ x>) I k 2 1 _ fc '2 : J_ > ) + m 2 X ') 2 



6 = 2A; ± fc^ (148) 
Then the functions in Eqs. fll43H145D are given 

Q 2 = a — b cos t 

5Q 2 = 5a, (149) 

where 

a = -{ai + a 2 ) 

5a = -(oi — a 2 ), (150) 

It is useful to display the matrix elements of the effective interaction in the form of 
tables. The matrix elements depend on the one hand on the momenta of the electron 
and positron, respectively, and on the other hand on their helicities before and after 
the interaction. The dependence on the helicities occur during the calculation of these 
functions E(x, kj_; X±, \2\x', k'j_; X[, X' 2 ) in part I and G(x, k±; Ai, \2\x r , k' ± ; X[, A' 2 ) in part 
II as different Kronecker deltas ||22|| . These functions are displayed in the form of helicity 
tables. We use the following notation for the elements of the tables 

Fi(l,2) -> E i {x,k ± )X',P ± )) Gi(x,k ± ;x',k' ± ) (151) 

Also we have used in both cases for the permutation of particle and anti-particle 

F;(x, k±; x', k' x ) = F 3 (l - x, -k±; 1 - x', -k' ± ) (152) 

one has the corresponding for the elements of arbitrary J z ; in the case when the function 
additionally depends on the component of the total angular momentum J z = n we have 
introduced 

Fi(n) = Fi(-n) (153) 



C.l The helicity table 

To calculate the matrix elements of the effective interaction in the exchange channel we 



use the matrix elements of the Dirac spinors listed in Table (|5|) p2 |. Also the following 
holds v X '{p)Tv x {q) = u x {q)l a u x >{p). 

We introduce for the matrix elements entering in the effective interaction Eqs. ( |143 - 

2E^(x,k ± ;X 1 ,X 2 \x',k' L ;X' 1 ,X' 2 ) = {^l u )9„, (154) 
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M 



7 



7" 



7l 



7l 



=u(k',\')Mu(k,X) 



26* 



k+k'~ 



(m 2 + k ± k' ± e +lX ^-^) ^ - mA (j^e^' - A; ± e + ^) 



I 1 



Table 5: Matrix elements of the Dirac spinors. 

with ( 7 ^)^ = |(7 + 7") + - (7?) - (7 2 2 ) and 

2E^(x,k ± ;Xi,X 2 \x',k' ± ;X' 1 ,X' 2 ) = (yy^^-L 

with ( 7 M 7 1/ )?7 M ?7,/ = (7 + 7 + ); where 

(w(x, A;_l; Ai) 7 ^m(x', A^)) (u(1 - x', — /c^; A' 2 ) j u v(l - x, -k±; A 2 )) 



<J(l-x)(l-x>) 



These functions are displayed in Table 



final : initial 


(Ai, A 2 ) =TT 


(Ai,A' 2 ) =U 


(A'i, A 2 ) =it 


(a;,a 2 ) =11 


(Ai,A 2 ) =TT 




£5(1,2) 


£3(1,2) 





(Ai,A 2 ) =n 


£5(2,1) 


£2(1,2) 


£4(1,2) 


-£3(2,1) 


(Ai,A 2 ) =|T 


£3(2,1) 


£4(1,2) 


£2(1,2) 


-£3(2,1) 


(Ai,A 2 ) =|| 





-£3(1,2) 


-£3(1,2) 


£1(1,2) 



(155) 



(156) 



Table 6: General helicity table defining the effective interaction in the exchange channel. 



The matrix elements El (1, 2) = E\ n> (x, k±; x', k' ± ) with n = 1 and n = 2 for 'g' and 



77^77^ terms, respectively, are the following 



E^ixM-x'^) 



m 



+ 



1 



k 1 fc'. 



-i(<p-tp ) 



xx' (1 — x)(l — x')y xx'(l — x)(l — x') 
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i?2 ^ (x, k±; x , k ± ) 



m 



xx' (1 — x)(l — x 1 ) 



k'l- 



'x(l — x) x'{l — x') 



+ k±k\ | — + 



xx 



(1 -x){l-x') 



E^ ^ (x, k±] x , 



m 

XX 
2 



l-x' , 



■ ~ 1 — x 

r\2 



up 



-m 



(x — x') 



xx'(l — x)(l — x') 



(157) 



and 



E{ '(x, k±;x' 

C.2 The helicity table for arbitrary J z . 



(x — x') 2 



(158) 



Following the description given in the main text Eq. ( p.20 ) we integrate out the an- 
gles in the effective interaction in the exchange channel. For the matrix elements of 
the effective interaction for an arbitrary J z = n with n 6 Z we introduce the func- 
tions G(x, k±; Ai, A 2 |x', k' ± ; \[, A' 2 ) = (x, k±; J z , Ai, A 2 |V r e fj|x / , k' ± ; J' z , \[, A' 2 ) in the exchange 
channel and obtain the helicity Table (0). 



final : initial 


(Ai,A' 2 ) =TT 


(Ai,A 2 ) =n 


(A'i, A 2 ) =|T 


(Ai,A' 2 ) =I| 


(Ai,A 2 ) =TT 


Gi(l,2) 


G*(l,2) 


G 8 (l,2) 





(Ai,A 2 ) =TI 


G$(2,l) 


G 2 (l,2) 


G 4 (l,2) 


-G 3 (2,l) 


(Ai,A 2 ) =|T 


G 3 (2,l) 


G 4 (l,2) 


G 2 (l,2) 


-G$(2,l) 


(Ai,A 2 ) =|I 





-G 3 (l,2) 


-G$(l,2) 


Gi(l,2) 



Table 7: Helicity table of the effective interaction for J z = ±n, x > x'. 
Here, the functions Gj(l, 2) = Gj(x, x', k' ± ) are given 



Gi(x, k±;x',k' ± ) 



i i 
m m 

+ 



+ 



xx' (1 — x)(l — x') 
k i k\ 



xx'(l — x)(l 



Int(\l - n\) 



25a 



aMM) ~ 7 K5 Ird(\l - n\) 

— X ) (x — x y 



G 2 {x, k±; x', k ± 



k'l 



a +-7TS+ TTT a Mini) 
') j x{l — x) x'{l — x') ) 



\ \xx' (1 — x)(l — x') 
+ fc^ - „|) 4- _ xl) lnt(\l + „|)) 
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—Int{\n\) 



(x — x'Y z 

m ( \, , .... , 1 — x' 



G 3 (x, k±; x 1 ', k ± ) = [ k\_Int{\l — n\) — k± Int(\n\) 

xx' \ 1 — x t 

G A {x,k ± -x\k^) = - m 2 — J^-^l— -^/nt(|n|) 



we define 



a f 27T cos nt 
I{n;a,b) — ———r / dt- 



2tt 2 Jo a — b cos t 

then in Eq. ( pL59| ) the following functions are introduced 
(1) Exponential cut-off 

Int(n) = I(n; a, b) 
Int(n) = 0, 



(2) Gaussian cut-off 



Int(n) = Kel(n; a + iSa, b) 
Int(n) = lml(n; a + iSa, b) 



(3) Sharp cut-off 



Int(n) = 8(—5a)I(n; a — 5a,b) + 8(Sa)I(n; a + Sa,b) 
Int(n) = 6(—5a)I(n; a — 5a, b) — 6(5a)I(n; a + 6a, b) 



also a + 5a = a\ anda — 5a = a 2 . 
Explicitly is used 



27r , cos nt 1 (a — \Ja 2 — b 2 \ 

dt = 2tt - 

a — bcost \Ja 2 — b 2 \ b J 

2n sin nt 

dt = , 

o a — b cos t 

where a can contain imaginary part as in the case of gaussian cut-off. 
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D The matrix elements in the annihilation channel 

We repeat the same calculations for the matrix elements of the effective interaction in the 
annihilation channel. Annihilation part of the effective interaction can be written 

V eS = -£i(YY)C»„ , (165) 



where one has 



V I M 2 + M 2 ) v+ 2 



a 



in the frame p± = 0.0 Explicitly the annihilation part of the effective interaction for 
different cut-offs Eq. fll3lf) - Eq. dl4l|) is given 
(1) Exponential cut-off 



= % , (167) 



(2) Gaussian cut-off 



r M 2 w 



(3) Sharp cut-off 



± j Q{Ml-Ml) 6{M 2 -M 2 a ) \ Wv 
^ = V I ^ + ^ J - (169) 



where p + = + is the total momentum; and (7^7^) for annihilation is defined in 
Eq. (|92"D. The functions present in Eq. ( |167| )- Eq. (|169| ) are given in the light-front frame 

M 2 = k '^ + m2 
x'{l-x') 

k 2 +m 2 , 

Ml = -j- r, 170 

x(l — x) 



2 Indeed (j^j^g^ = |(7 + 7") + |(7^7 + ) + (t'VOs^; therefore it holds 

_l -Py) +Vv(P» -vj) _ P^_ 

9^v — 9uv + + ,2 r l^ 1 lf > 

P p + 

The 4- momentum of the photon p^ in the ^-channel can be written p^ = p lfl + p 2fl — ri^D a /2 = p\^ + 
P2fj, — VnDb/Z with D a ,Dt defined in Eq. (pOj). The Dirac equation (pi + P2) ( uu(pi)7 Al, v(p2) = allows 
then to write p^u{pi, Xi)j ,J 'v(p2, X2) — —M^l 2p + )i]^u(pi, Xi)j tJl v(p2, A2). Thus, when p± = 0, one has 

-> — ~T M i 

where the arrow means that this tensor should be contracted with (7^7^) in the annihilation channel. 
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we remind also 

M 2 = \{Ml + Ml) 

SM 2 = \{Ml-Ml) (171) 

Note that the energy denominators of the effective interaction in the annihilation 
channel do not depend on the angles ip, if'. 



M 


—L=v(k',X')Mu(k,X) 


7 + 




T 




6i x , - mX ( 


k' ± e +iX<p ' + k ± e +iXlp ) 4] 


7l 


( k '± -<V , k ± +iX>p^ 
[k'+ k+ J 


6 X _ X , - mX ( 


,k'+ + k+) 




7i 


iX (^ e -V_^ e+ a^A v _. m 


{l7+ + k+, 





Table 8: Matrix elements of the Dirac spinors. 



D.l The helicity table 

For the calculation of matrix elements of effective interaction in the annihilation channel 
we use the matrix elements of the Dirac spinors listed in Table @ [22] . Also the following 
holds (v X '(ph a u x (q)) + = u x {q)^ a v x >{p). 
We introduce 

2HM(x, k ±] Ai, X 2 \x', k' ± ; X[, X' 2 ) = (WOfli = -(if) - ( 7a 2 ) 
2H^(x,k ± ;XuX2\x',k' ± ;X / 1 ,X / 2 ) = (172) 

where 



(w(l-x / ,-^;A^)7^u(a; / ,A;l;A' 1 ))(M(x,A; ± ;Ai)7^(l-a;, -fcj.; A 2 ) n - Q , 

\1 1 ) = / . / . ( 173 ) 

Jx'(l — x') Jx(l — x) 
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final: initial 


(Ai,A' 2 ) =TT 


(a;,a 2 ) =n 


(A'i, A 2 ) =|T 


(A'!,A 2 ) =|I 


(Ai,A 2 ) =TT 


#i(l,2) 


H 3 (2,l) 


# 3 *(2,1) 





(Ai,A 2 ) =n 


H 3 {1,2) 


H*(l,2) 


ff 4 (2,l) 





(Ai,A 2 ) =IT 


#|(1,2) 


Hi(l,2) 


#2(1,2) 





(Ai,A 2 ) =I| 















Table 9: General helicity table defining the effective interaction in the annihilation chan- 
nel. 



These functions are displayed in the Table @. 

Here, the matrix elements H^ n \l, 2) = H^ n \x, k±; x', k' ± ) are the following 



and 



%aj _L iAj / \ »X/ -L. \Xj 

e i(<p-(p) 



H^\x, k±; x', k ± ) = —m\i(— 



x 1 — xj 1 — x' 



e i(<P~P ) 



H^ix^r^x',^) = k ± k' ± I _ 1 (174) 



H{ 2) (x, fc±; x', jf ± ) = flf fc±; a/, k' ± ) = 

i/" 2 (2) (x, jfej.; x', k' ± ) = Hf\x, k ± ; x', k' ± ) = 2 (175) 



D.2 The helicity table for | J z \ < 1 

The matrix elements of the effective interaction for J z > F(x, k±; Aj, A 2 |x', A' 1; A 2 ) = 
(a;, J z , Ai, A 2 |K//|a ;/ , ^jj A' 1; A 2 ) in the annihilation channel (the sum of the gener- 
ated interaction for J z = +1 and instantaneous graph for J z = 0) are given in Table ([TOD - 

The function F;(l,2) = F, the following; 



J-? / , / ? ' \ "1 ^ 2 x 

F ^> k ^ x > k J = nnxx>(l-x)(l-x>) 5 ^ 



F 2 (x,k ± -x',k' ± ) = "(1^5^ + 25^) 

7T \\l XX J 

F 3 (x,kr,x,k ± ) = -^ xl[1 _ xl)l _J w 
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final: initial 


(Ai, A' 2 ) =TT 


(a;,a 2 ) =n 


(a;,a 2 ) =it 


(Ai,A' 2 ) =|| 


(Ai,A 2 ) =tt 






F 3 *(2,l) 





(Ai,A 2 ) =n 




F*(l,2) 


^4(2,1) 





(Ai,A 2 ) =|T 


F 3 *(l,2) 


^4(1,2) 


^2(1,2) 





(Ai,A 2 ) =|| 















Table 10: Helicity table of the effective interaction in the annihilation channel for J z > 0. 



^ i i i i' \ a ( I k ± k' ± . n . N 

F 4 (x,k ± ;x,k ± ) = -( v -^^T^)^U + 2 ^ 



(176) 



where we have introduced 
(1) Exponential cut-off 



(2) Gaussian cut-off 



(3) Sharp cut-off 



1 



l l 



1 _ Ml + Ml 
Q Af 4 + Mt ' 



O(Ml-Ml) , ^(M 6 2 -M a 2 ) 



(177) 



(178) 



M„ 2 



+ 



M 6 2 



The table for J 2 = — 1 is obtained by inverting all helicities, i.e. 

F{ J z = +1; Ai, A 2 ) = -X 1 F(J Z = -l;-X 1 , -A 2 ) , 



(179) 



(180) 



The matrix elements of the effective interaction in the annihilation channel are nonzero 
only for \J Z \ < 1 due to the restriction on the angular momentum of the photon. 
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